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ABSTRACT 


A technique to remotely determine the attenuation 
coefficient (a) of the water was investigated. The back- 
scatter energy (0 = 180°) of a pulse laser (X = 440 -• 

660 nm) was found directly related to the water turbidity 
for waters characterized by a <_ 6m ^ . The greatest sensi- 
tivity was found to exist at 440 nm. For waters whose 
turbidity was adjusted using Chesapeake Bay sediment, the 
sensitivity in determining a at 440 nm was found to be 
approximately 5 — 10%. A correlation was also found to 
exist between the water depth (time) at which the peak 
backscatter occurs and a. 


ACKNOWLKl>ftf:MKN'J .• 


The authors world like to acknowledge the* able 
assistance of technician, Mr. Al Boshoff, who helped 
in performing the experiments in SPARCOM'S environmen- 
tal tank. 

Dr. Peter Mumola, NASA Langley, and Mr. Les Dunn, 
EPA, Las Vegas are gratefully acknowledged for serving 
as project monitors. 


T 


I. Introduction 

Lasers 1 ' 2 ' 3 have been used to remotely measure proper- 
ties of the water. Precise interpretation of the return 
laser signal (s) requires knowledge of the absorption and 
scattering coefficients of the water as a function of wave- 
length. 

The quantity and type of sediment suspended in water 
determines the turbidity of the water and the transmission 
of electromagnetic energy in water. The results of a pre- 
vious study showed that there was a direct correlation be- 
tween the quantity, type of sediment and the optical attenua- 
tion coefficient (cx) of water. The object of this study 
was to determine the feasibility of using the laser back- 
scatter signal to remotely determine the attenuation co- 
efficient of 'the water. If proven feasible, a device based 
on this measurement could be an integral portion of an air- 
borne laser bathymetric surveillance system. _In addition; 
it could be used to map the discharge of heavy "sediment laden 
streams, etc, into larger bodies of water. 

II. Optical Properties of Water 

The attenuation coefficient (a) of the water is given 
by the equation 

a (X) = a (X) + s (1) 
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civ dif fcrentinl scattering volume 

dl(0) « radiant intensity at an angle 0 

0 angle with respect to forward direction of 

incident beam. — 

The scattering coefficient (s) is obtained by integrating 
the volume scattering function over the solid angled ft 
(Equation 4) . 

ir 

s = Jo (0) dft = 2 nr J o(G) sinO dG (4) 
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III. Experimental Measurements 

The laser backscattered signal was measured in Spar- 
corn's Environmental Water Tank using two sediments having 
different values of a/s*. The two sediments that were used 
to control the turbidity in the water tank were quartz (a/s - 
0.06) and Chesapeake Bay sediment (a/s = 0.20) l! The wave- 
length of the laser light was varied from 440 to 660 nm in 
intervals of 20 nm. Figure 1 is a schematic diagram of the 
experimental set-up that was used to measure the backscattered 
signal . 

An additional set of experiments were completed to cor- 
relate the attenuation coefficient with Secchi disk measure- 
ments . The latter are commonly used for measurement of 
turbidity. Appendix A contains some details on the use and 
limitations of Secchi disk measurements. 

The beam from the pulsed dye laser entered the Environ— 

*a/s ratios were measured at 540 nm. 

(3) 
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Figure Is A schematic diagram of the experimental set-up that was used to 
measure the backscattered signal. 





mental Tank from above the water a:: .shown in Figure 1, The 
phot omul kiplior tube (PMT) viewed the laser baekscattered 
siqnal from 1.7 meters above the. water . Detection of the di- 
rect surface reflection was avoided by tiltxng the detector 
approximately 2° with respect to normal incidence of the la- 
ser beam. 

I he effect of possible interfering surface and bottom re- 
flections on the backscatter signal was investigated. This 
was accomplished by placing highly reflective surfaces both 
on the water surface and the bottom. It was discovered that 
the -slight tilt on the detectors, as previously mentioned 
was sufficient to avoid reflections from the highly specular 
water surface. For clear shallow water (a ■ <_ 0.5 m 1 ) the 
backscatter is small in relation to the bottom signal. How- 
ever, for values of a >_1.0 m" 1 the backscatter becomes the 
predominate term and is easily distinguished from tua bot- 
tom signal. 

The present set of experiments was performed for quies- 
cent waters. Future experiments should be designed to in- 
vestigate the effect of capillary waves on the backscatter 
signal . 

A series of experiments were performed to measure the 
shape, the maximum amplitude and the depth of water from 
which the backscattered laser pulse is generated. These 
measurements were also performed as a function of the de- 





( 5 ) 


tcctor's fic*] cl of view. Ac previously mentioned, the wave- 
length of the laser war. varied from 4 4 0 - GGO nm. The re- 
sults of the quartz measurements are summarized in Tables I 
and II for an input laser power of 1 watt. The results of 
similar type experiments using Chesapeake Bay sediment are 
summarized in Tables III and IV. Figure 2 shows sampling 
backscattered signals for a wavelength of 540 nm from tur- 
bidities characterized by a values ranging from 1.0 - 6.0 m~ 1 . 
These results were obtained for quartz induced turbidities 
with a 200 mrad. field of view detector. However, the shape 
of the curves are very similar to those obtained for the 
Chesapeake Bay sediments. They were found to be independent 

of the detector field of view. The amplitude of the pulse_ 

shown in Figure 2 cannot be compared as they were observed 
using various detector sensitivities. 

The relation between the amplitude of the peak backscat- 
tered signal and the attenuation coefficient a is demon- 
strated in Figure 3 for the case where Chesapeake Bay sedi- 
ment is used to control the water turbidity. Data are shown 
for four different wavelengths (X = 440, 480, 500 and 560 nm) 
for a detector having a 200 mrad. field of view. Figure 4 
shows similar type data for the case where the turbidity of 
the water was controlled by quartz particles. The curves 
shown in Figures 3 and 4 represent a best fit to the experi- 
mental numbers. It should be noted that in all cases the a 
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and a valuer; have been measured at 5 4 0 run. .... M o n : ; u r e i ne n L 
are- currently in progress to determine the wavelength do- 
pondoncu of both a and a in the wavelength region of 440 
to G60 nm. These results will bo used to correct the u 
values given, in this report. The shape of the curves in 
Figures 3 and 4 are basically the same except that the 
slopes are appreciably steeper for the natural sediment, 
thereby resulting in better sensitivity .- 

In order to compare the behavior of the various curves 
shown in Figures 3 and 4 a quantity R^ has been defined as 

peak backscattered signal @ a = 6.0 m ^ 

R = Ij (5) 

x peak backscattered signal @ a = 0.5 n 

where x represents the particular particles used to simulate 
the water turbidity, i.e., 

x = s (Chesapeake Bay sediment) 

= q (quartz particles).... 

Table V gives the values of R (X) and R (X) for wave- 

s q 

i 1 lengths of 440 to 660 nm. 


X (nm) 

440 

480 

500 

560 

580 600 

620 

660 

to 

10.0 

5.00 

5.44 

4.86 

6.00 

- 

yx> 

3. 89 

3.13 

2.64 

2.74 

3.80 3.50 

3.50 

2.40 


Table V: Ratio of peak backscattered signals for differ- 

l ent wavelengths. 


It ^ 
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RELATIVE AMPLITUDE OF PEAK BACKSCATTERED SIGNAL 



Mm" 1 ) 


Figure 3: The amplitude of the peak backscattered 

signal is shown as a function of a for 
various wavelengths. The turbidity was 
adjusted using Chesapeake Bay sediment. 



Figure 4 : The amplitude of the peak background 

signal is shown as a function of a 
for various wavelengths. The turbidity 
.was adjusted using quartz. 
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The* results given in Table V .chow 
a curve exhibit s a posit: j ve slope 
This means that the backseatter si 


that the; backseatter/ 
for all wavelengths, 
gnal from any wavelength 


in the 440 - 6C0 nm region could bo used as a measurement 
of a. For both quarts . and .natural sediment the slopes of 
the curves are approximately the same for all wavelengths, 
except at 440 nm. The slope of the backseatter curve at 
the latter wavelength is substantially greater. A laser 


operating at this wavelength would be the most sensitive 
to changes in a, thereby making it the best a detector. 

The sensitivity o i this technique, using the wavelength of 
440 nm, is estimated from Figure 3 to be approximately 5 - 
10%. A laser operating at still a lower wavelength may 
result in increased sensitivity. 

It has been established that in some underwater appli- 
cations the use of polarization techniques reduces the back- 
scattered signal and improves the contrast. 5 An experiment 
was performed to measure the effect of polarization on the 
backscattered signal as a function of turbidity. The laser 
beam was linearly polarized using one polarizer in front of 
the laser output while a second polarizer was placed in front 
of the detector. By changing the relative angle of the two 
polarizers, the amplitude of the backscattered signal could 
be maximized or minimized. It was determined expcrimently 
that there was no correlation between the ratio L 


i 
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amplitude of the maximum signal 
amplitude of the minimum signal 
and the turbidity of the water. 

IV. Analytical Model 

The shape of the laser backscattored signal and the 
time at which its peak response occurs in the water is 
predicted using a one dimensional model for the experi- 
mental geometry shown in Figure 5. This model is designed 
to examine the primary featiares of the backscattering prob- 
lem for the case of a pulsed laser. Additional calcula- 
tions are needed to develop a three dimensional analytical 
model for the prediction of the detailed shape and the mag- 
nitude of the backscattered signal. 

At a time t - 0, a short pulse is propagated towards 
the water and impinges on the air/water interface at a - t i me , 
t = D/c where D is the distance- between the laser/detector 
system and the air/water interface. Consider the case for 
which the front edge of the pulse is at position x and the 
entire pulse is in the water. The parameter t fixes a point 
within the pulse. The amplitude of the signal detected at 
a time 2t,f(2t), neglecting the difference between the velo- 
city of light in air and water, is given by 

L/2c 

f ( 2 1 ) = J g(2x) P (X ; t-T ) dx (6) 

T = 0 
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t > (D+U/c 


anu 


lf/2c 


f(2t) = 


f g( 2t) P (X;t-t) dt 


(7) 


J 

x=0 


for 


(D+L)/c > t > D/c and 0 < L < L 


where 


L is the length of the pulse, 


f (2t) = 0 for t < D/c 


( 8 ) - 


g( T ) is the power distribution function for the pulse. 
The peak amplitude of the pulse occurs at V such that 
g ( x ) si and the peak power is ^ ^ 

the response function of the sediment which is repo 
for the generation of the backscattered signal. 

The response function P(X;x) is derived for a line ele 
ment located at a point x as shown in Figure D 


dP(x.-x) = (p./D e- 2 “> ( *- D) o(X,.) aa 


(9) 


which upon integration yields 


- (P./D e- 2 “x. (x - D) o(x,.J 2 " (1 - 7^T7T7I 


) (1 


P ( X ; x) 


whor „ a i, the radius of the detector. The exponential 
tone account.'; for the attenuation of the nic;n. I rn the wn 
ter while a(X;ir) is the volume scattering function at 1«0 

degrees « 

Inserting equation (10) into equation (6) results in 
thr . following expression for the detected signal « - 
function of a, and wavelength, i.e., 


f(2t) » (Po/L> 2lTC o(X?1t) e 


-2<x*e(t--D/c) I(t) 


where 


lit) •- 


A/2? (t-D/c)/2 / 

f d , giw A CT i- 

-0 1 


C (t-T) 


(a 2 +c 2 (t-r ) 2 ) 1/2 j 


with y * L/C. The upper limits of I (tj * efer to the 
(y/2 , when the whole pulse is in the water or when part of 
the pulse (tt-D/oI/J) i- m the water. The ...odified Gaus- 
sian pulse shape given by equation (13), is assumed for the 
power distribution function ,(t> which is inserted into 

equation (12) . 


g(x) = 


(2T/Y-U - e" 1 ) 


-it u 


The return laser signal has been calculated by »u»orioally 
integrating equation (12) . Appendix 1) describes the prograi 


I ^ 


— *m V &1: 
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for the manor leal integration. The ror.u.l 1 :: of the calcula- 
tions for D 1.7 motor.'? , a detector raid in:? (a) of 0.5 era 
and o. vain o.‘» ranging from 1 - 4 m ^ arc shown in Figuro 6. 
The value of o(A;r) was assumed to ho equal to 1 since the 
value of this function for the simulated waters was unknown. 
As the value of a is increased the peak of the backscattered 
signal shifts toward the surface of the water. It should be 
noted .that, these curves only show the shift of the peak with 
time. The relative amplitudes have not been normalized and 
therefore the various curves cannot be compared in amplitude 
Table VI lists— the time displacement of the pulse with res- 
pect to the water surface 


a (m -1 ) 

1 

2 

3. 

4 

Time displacement 
of pulse (nsec) 

5.8 

4.7 

3.9 

3.4 


Table VI: Time displacement (nsec) of the backscattered 

pulse for various values of a and y = 6.0 nsec. 

The experimental pulse shapes shown in Figure 2 show 
the same— time_shif t as those displayed in Figure 6. Fig- 
ure 7 shov;s the comparison between the experimental and 
calculated time shift as a function of a. Two calculated 
curves are presented; y = 6 nsec, y = 5.5’ nsec. One ob- 
serves that a definite correlation exists between the wa- 
ter depth at which the peak backscattcring occurs (time 



i 


I 





Figure 7 : Comparison between the experi- 

mental and calculated time dis 
placement of the peak back- 
scattered signal. 
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displacement: ) and the wator turbidity (a). In .addition, 
good agreement: is noon to cxa.nl befwt on (lie calculated and 
exper.i mental points . 

V. Summary 

A series of experiments were performed in which the 
laser backseatter energy was measured as a function wave- 
length and water turbidity. A tunable dye laser was used 
to vary the laser emission from 440 to 660 nm, The turbid- 
ity of the water was changed in discrete steps for values 
of the optical attenuation coefficient a < 6 This 

was obtained by adding quarts and Chesapeake Bay sediment 
to the water in SPARCOM'S environmental -water tank facility. 

The results showed a definite dependence of the back- 
scatter energy with the attenuation coefficient (a) at all. 
wavelengths. However, the strongest dependence of the back- 
scattered signal oh a was observed at 440 nm. The depend- 
ence on a was the same for both sediments, however, it was 
more pronounced for the Chesapeake Bay sediment. At an 
operating laser wavelength of 440 nm it is estimated that 
the attenuation coefficient (a) can be determined to within 
5 - 10% in Chesapeake Bay type waters. 

It was also shown that the peak of the backscattered 
energy occurred at a distcnce in the water corresponding to 
a two-way time path of 2.5 - 5.5 nsecs. The peak shifted 
towards shorter times, i.e. , shallower depths as the water 
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turbidit y increased . A ono di mansion.'! 1 analytical model 
was formulated which war; able to predict belli the absolute 
time peaking and shift of the bacdcsealtcri d signal. 

The initial calculations and measurements have shown 
that it may be feasible to use the laser signal backseat- 
terec'i from the near surface of the water as a remote sensing 
tool of -water turbidity. 

Although the backscatter is very sensitive to a, it may 
not be possible to design a system to measure -he absolute 
value of a by detecting the laser backscatter at a single 
wavelength. However, it appears that a technique could be 
developed— that is based on detecting the ratio of the laser 
backscatter at two or more wavelengths. Where applicable, 
a ratio technique can be extremely useful in that the major- 
ity of the system parameters cancel out. Additional measure 
ments at various wavelengths, along with more detailed cal- 
culations are required in order to arrive at a ratio tech- 
nique that can be used to determine a on an absolute basis. 
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Appendix A - Measurement of Water Turbidity 

The optica] properties of turbid water can be uniquely 
determined by urc of tv.'o active 1 i g h t measuring instruments; 
an alpha (a) meter and an absorption (a) meter. The cx-meter 
measures the total transmission loss of a light beam in v’ator 
due to absorption and scattering, while the a-metcr measures 
the loss due only to absorption. 

Instead of using the above devices, which are relatively 
expensive and tedious to use, oceanographers for years have 
relied mainly upon Secchi disk readings for determination 
of visibility. While the Secchi disk is inexpensive ,. it is 
also less accurate than the active devices and is much more 
susceptible to errors produced by weather conditions and 
operator interpretations. 

The- Secchi disk 6 is a diffuse white target of radius r 
and reflectance p D that disappear s-from view at a depth D 
in the water. The visibility of the disk is determined by 
the contrast between the disk and the background given by 

_ , , Disk Brightness - Background A 1 

Contrast =. - - - — — — * A. 1 

Background 

The background light -reaching the eye is the result of 
surface reflection and upwelling light. The background 
irradinnee from the surface (B ) is given by the equation 

5 
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where 


o = surface reflectance 
' s 


E 0 = incident intensity of light 

u = fraction of incident radiation scattered upward 
The background irradiance produced at a depth D is given by 
the equation 

“kD 


B d = u E 0 e ~ ( 1 -P s )/tt 


A. 3 


where 


k = diffuse extinction coefficient. 


The flux from the disk at depth D is 

— kn 

<t» D - Eo (1-P S ) p d d A 


A. 4 


where 


P D = reflection from the target at depth D 


and is attenuated by a factor of 

e” aD (l-p s )/(D 2 +r ? ) . 

Thus, the. total ..flux from the disk that reaches the eye of 
the observer .is given by 

v - (k+cx) D 

A. 5 


0 = 
D 


E 0 (1-P S ) 2 P D e 


it (D 2 +r 2 ) 
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From the equation A.l, the contrast at. 
given by 


C D = 


4> -B 
D D 


B 


D 


depth 1) (C D ) is 


A. G 


or 


C 


D 


P D “ U 

u 


A.? 


However, the contrast seen by the eye is given by 


eye 


= C D { 


B D , -ccD 

^ e 


A. 8 


.in which the bracket represents the attenuation of the 
contrast of a submerged object. . 

If the contrast threshold is chosen. to be 0.003, the 
sum of the attenuation and diffuse extinction coefficients 
as obtained from equations A. 2. A. 3, A. 6 and A. 8 is 


k+cx = 1/D In [333 (p D ~u) / (u+w) ] 


where 


w « p s tt/(1-p s ) 


A. 10 


When a glass bottomed viewer is used (p g = 0 ) , equation 
A. 9 becomes 


k+t = 1/D In 1333 (p d /u-1)] 


A. 11 


i ^ 


\C 
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For „ constant value of a, Holmes’ ha, determined that 


D (k+a) = constant ^9.4 


A. 12 


In addition, the experimental data has boon fitted to 
produce the equation 

a = In [100/(5. 3D-3. 4)] 


which is used to establish a relationship between k and a. 

The a-meter provides accurate and repeatable measurements 
for a variety of operating conditions. A light source is 
focused to form a diffraction limited beam that is propa- 
gated through 1 meter of water. The light is detected by 
a small aperature photodetector. However, the finite spread 
of the beam in the conventional a-meter produces a value of 
a that is 25% lower than that obtained using an. extremely 
narrow beam laser transmissometer . Another problem inherent 
with the conventional a-meter is the broad spectral 
ing which is used on the detector. For the present measure- 
ment a narrow band (10 nm) interference filter was placed 
in front of the photodetector to enable measurements to be 
made of a as a function of wavelength. 

The a-meter uses a ratio technique to measure the absorp- 
tion by the water. A point source of light is placed between 
two detectors having wide fields of view. Each detector is 
a different distance from the point source so that the dctectc 
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radiation of each detector is given by the equation 


P (i) ~ P e** aD i A. 14 

D s 

where 

P - power of source 

= area of detector — 

D i = distance between i t ^ detector and the point source 
i - 1 , 2 detector number 

a - absorption coefficient — . — 

The absorption coefficient is then determined from equation 
A. 15 

(D x -D 2 ) a = Jin { P D (1) /P D (2) } -2 In (D 2 /D 1 ) ] A.15 

This is accomplished electronically by a logarithmic ampli- 
fier. The values of and D 2 were chosen to be 1/3 and 
2/3 meters_respectively . Honey and Sorensen 8 have shown that 
accuracies better than 1% are achievable under moderate 
scattering conditions and accuracies of better than 7.5% are 
achievable for severe scattering conditions. 

A set of measurements were made to compare the perform- 
ance of the a-meter and the Secchi disk. This was accom- 
plished by using the a-meter to measure the attenuation co- 
efficient at 540 nm while observing a one foot diameter half- 
black half-white Secchi disk. The disk was moved away from 
the observer until it disappeared thereby defining the dis- 
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tance D. For small and moderate u'c this was done by moving 
the dick along the long dimen cion of the tank. For larger 
values of a the measurement: v/ac made by lowering the dick 
vertically into the water. In both cases a largo fluores- 
cent fixture was used to illuminate the disk. Shown in Fin 
ure A. 1 is a plot of the experimental__data along with Holmes' 
equation (A. 13). 

An inspection of this figure shows that for values of 

-1 

a £ 2 m , the difference between the measured and predicted 
values is Aot - 0.5 m For larger values of a this differ- 
ence becomes greater, i.e—, Act = 1.0 m” 1 . Possible sources 
of error are 1) the sun was the light source used in Holmes' 
experiments, whereas a standard fluorescent fixture was used 
in the present experiments, 2) differences in the intensity 
of the various sources and their spectral distributions. 


Appendix D - Computer Program 


The integrals expressed in equations (10) and (11) 
were integrated numerically using the Gnuss-Dimpson Rule. 
The integral 



dx 


can be represented by a summation 



dx = 


h 


n/2-1 


z 

m=r A 


| f (a+hki) + f (a+hk 2 ) } 


(C.l) 


for the values of n even and h=(b-a)/n where 


ki » l~l/V'3+2m k 2 = l+l//3+2m 


(C 


,. 2 ) 


The program used to evaluate equation C.l is given be- 
low. The parameters used in the calculation are explained 
after the computer program. 


Main Program 


C=. 29 97 
C2=C*C 
D=l. (.50 
A= .005 
A2=A*A 
ZIGMA=1 . 
P0=1. 
EL=1 . 5 


“'."H IPS- 
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GAMA- 6.0 
ALPHA 1 ' 1. 

N-10 

SN-N 

M=100 

HM-M 

CC= 2 . *3 . 14 16 GMA * V 0/E L 

DO 1 J=l,6 

WRITE (6,41) ALPHA, GAMA 
4 1 FORMAT ('0',' 7vLPIIA= ' , F 5 . 2 , 3X , ' GAMA- 
T=D/C 

DO 2 1=1, N 
S=0. 

AL= (T-D/C) -GAMA 
IF (AL) 111,112,113 

111 TAMAX=0.5* (T-D/C) 

GO TO 114 

112 TAMAX=0.5* (T-D/C) 

GO TO 114 

113 TAMAX-GAMA 
114' HTA=TAMAX/HM 

Al=-1. 577350 
100 TA=Al*HTA+2 . *HTA 
LL=M/2 

DO 30 JL=1 ,LL 
Fl=C2* ( (T-TA) * * 2 ) 

F1=A2+F1 
Fl=SQRT (FI) 

F=C* (T-TA) 

Fl=l . -F/Fl 

F2=EXP (2 . *ALPHA*C*TA) 

Gl=4 . *TA/GAMA-1 . 

Gl=-G1*G1 

G1=EXP(G1) 

Q=GXP (1 . ) 

QV=1./Q 

FM=Q/(Q-1.) 

Gl-FM* (Gl-QV) 

G1=G1*G1 

G2=CC*EXP (-2 . *ALPHA*C* (T-D/C) ) 

ENT=C*F1*F2*G1*G2 

S=S+ENT 

30 TA-TA+2. *HTA 

IF (Al.NE. -1.577350) GO TO 44 
Al=Al+2 „ /SQRT (3. ) 

GO TO 100 
44 S=HTA*S 


\ 
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,E13. 3) 


WRITE (6,-12) T , S 

4 2 FORMAT ( ' , ' T- ' , E 1 1 . 3 , 2 X , ' P>,~ ' , El 3 . 3 ) 
2 T=T3 GAMA/HN 
1 ALPHA- ALPHA-11 . 

CALL EXIT 
END 


Program Parameters 

C = velocity of„_light = 0.2997 m/nsec 

D = distance between detector and water surface in 
meters 

A = radius of the detector in meters 
EL = pulse length in meters 
P0 = total power in watts 

ZIGMA = volume scattering function in the backwards 
direction 

GAMA = pulse duration in nsec — 

ALPHA = attenuation coefficient of water in inverse 
meters 

N = number of intervals that the pulse has been 
divided 

M = number of terms in the summation 
TA = integration variable 

TAMAX = upper limit of integration 


ENT = integrand 
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